Different intermittency for longitudinal and transversal turbulent 

fluctuations 



Siegfried Grossmann, Detlef Lohse, and Achim Reeh 
Fachbereich Physik der Universitdt Marburg, 
Renthof 6, D-35032 Marburg, Germany 
(February 5, 2008) 

Scaling exponents of the longitudinal and transversal velocity struc- 
ture functions in numerical Navier-Stokes turbulence simulations with Taylor- 
Reynolds numbers up to Re\ = 110 are determined by the extended self sim- 
ilarity method. We find significant differences in the degree of intermittency: 
For the sixth moments the scaling corrections to the classical Kolmogorov 
expectations are 5£q = —0.21 ± 0.01 and 5£q = —0.43 ± 0.01, respectively, 
independent of Re\. Also the generalized extended self similarity exponents 
Pp,q = fi£, p /8£, q differ significantly for the longitudinal and transversal struc- 
ture functions. Within the She-Leveque model this means that longitudinal 
and transversal fluctuations obey different types of hierarchies of the mo- 
ments. Moreover, the She-Leveque model hierarchy parameters (3 L and (3 T 
show small but significant dependences on the order of the moment. 
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I. INTRODUCTION 



One of the central issues in turbulence theory has always been whether the velocity 
structure functions deviate from Kolmogorov's classical expectation [0,0. For many years 
the community focused on the longitudinal structure function 

D L p (r) = ([(u(x + r) - u{x)) ■ e L r f) = ({v L {r)f) (1) 

whose inertial subrange (ISR) scaling exponents we define as £p . (Here, is the unit 
vector in r direction.) The reason for this was that in many experiments Taylor's frozen 
flow hypothesis [l|-|3| had to be employed and therefore, the transversal structure functions 

Dj(r) = ([(«(« + r) - „(*)) ■ e T r f) = {{v T {r)f), (2) 

being a unit vector perpendicular to r, could not be obtained. We denote the scaling 
exponents of -Dj(r) as (p. 

A priori, there is no reason to expect = (J for general p. The PDF of the longitu- 
dinal velocity difference v L {r) is skewed because information from x to x + r is conveyed 
by the velocity difference itself, and odd moments of v L (r) thus do not vanish, Dp ^ 0. 
In particular, for homogeneous and isotropic turbulence the third order longitudinal struc- 
ture function is connected to the second order one by the Howard-v. Karman-Kolmogorov 
structure equation 

Di(r) = ~er + 6u^ : D^r). (3) 

Here, e is the mean energy dissipation rate and v is the kinematic viscosity. The PDF of 
the transversal velocity difference v T (r), on the other hand, is symmetric and consequently 
all odd order moments vanish, -Dj(r) = 0, p odd. Only the second order longitudinal 
and transversal structure functions are expected to scale the same way |l],f|,|5| in the large 
Reynolds number limit since for isotropic flow incompressibility implies 0] 

Dl{r) = D L 2 {r) + r -±-D L 2 {r). (4) 

Recently, multi-probe and optical techniques made the transversal structure functions 
experimentally accessible @,|6|-10|. In addition, also recent numerical simulations of decaying 



turbulence |Tl| have focused on the difference in scaling of D^(r) and Dj(r). 

The results of all this work seem to be contradictory. While all authors agree that (i) 
both (p and Cj show significant deviations from the classical expectation p/3, and that (ii) 
(p is well fitted by the She-Leveque model ("SL model") fL2"f , saying that 

C P = f-^o(f(l-/5 3 )-(l-/5 p )) (5) 

with C = 2 and (3 = (2/3) 1/3 , they disagree on whether = (J or ^ (J- 

Van der Water's group |7],|8[ finds that in shear flow the (moduli of the) transversal 
intermittency corrections 5(p = £j — p/3 are significantly larger than the longitudinal ones 
5(p = (p — p/3, p > 3. The same is found for jet flow turbulence by Camussi and Benzi [10 
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and for decaying (numerical) turbulence by Boratav and Pelz [11 1. On the other hand, the 
recent experiments by Camussi et al. || and Kahalerras et al. [jl3| did not give significant 
deviations between S( p and S(p and experiments by Noullez et al. |§ found 5C^ comparable 
to the 5(p found in other experiments. For a simple quantification of the intermittency 
corrections it is common to give 5( 6 . We do so in table | for the addressed experiments and 
simulations. 



ref. 


flow 


Rex 


"^6 




remark 


van der Water et al. []|] 
Camussi and Benzi [10 

Camussi et al. || 
Boratav and Pelz [|ll]] 
Noullez et al f§ 
this work 
this work 


shear flow 
jet flow 
wind tunnel flow 
decaying numerical flow 
jet flow 
forced numerical flow 
forced numerical flow 


up to 600 

250 

37 
~ 100 
up to 600 

110 

70 


0.18 - 0.20 

0.25 
0.19 ±0.03 

0.23 

0.21 ±0.01 
0.22 ±0.01 


0.27-0.31 

0.38 
0.22 ±0.03 

0.43 
0.25 ±0.10 
0.43 ±0.01 
0.43 ± 0.01 


normalized 
ESS 
ESS 
ESS 
ESS 
ESS 
ESS 



TABLE I. Intermittency corrections <5£q and (see text for the definition) for the longitudinal 
and transversal sixth order structure functions for various numerical and experimental flows. 



We would like to caution the reader of a too simplistic interpretation of table 1. The 
detailed definitions of the #Cp' T slightly differ from experiment to experiment. In all low 
Reynolds number experiments ||, PHTTH the scaling exponents could only be determined by 
employing the extended self similarity (ESS) method introduced by Benzi et al. [fL4| . In 
this method the structure functions D p (r) are plotted against D%(r), where D%(r) is the 
third order structure function defined with the modulus of the velocity difference and is 
experimentally found to scale with roughly the same exponent as D 3 (r), which includes 
the sign of v(r) [pT5|| . The ISR scaling exponent of such a plot is henceforth denoted as 
£p = Cp/Q- As in the ISR = 1 according to the structure equation ([3|), one expects 
£p = Cpi if Q L = C^- In experiments, however, there are always small deviations from 
Q L = 1; therefore, in principle and could slightly differ, see also ref. []TB|] . For the 
transversal structure function D^{r) there is no known relation as equation (FJ). Indeed, 
van der Water's group finds Q T = 1.08 |J and they give normalized scaling exponents Cp/Cs 
which we also call £ p . To date, there is no strict theoretical argument why ESS works so 
well. Note that for the present numerical simulation it does not work for odd order structure 



functions, taken without the modulus [|16[ . 

What is the origin of the differences between the results reported in table 1? One may 
think that it is the different geometry of the flows which causes the differences, in particular, 
a different strength of the local shear and of the anisotropy in the flow. Indeed, the shear in 
the flow of ref. [H] is considerable and it is known that shear destroys ESS |TT|,|T8|. On the 



other hand, at least <5£ 6 was found to be remarkably independent of different flow geometries 
p~9|] . Also, the numerical flow of ref. [Tl[] which clearly shows <5£g ^ is highly isotropic. 



But it is decaying which in ref. [20| is speculated to be a possible origin of the observed [11 



discrepancy between longitudinal and transversal intermittency corrections. 

In this paper we set out to determine the scaling exponents of the longitudinal and 
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transversal structure functions for forced, statistically stationary numerical Navier-Stokes 
turbulence up to Re\ = 110. Our motivation is to contribute to clarify the contradictionary 
picture reflected in table 1. It is of general interest for the understanding of intermittency 
whether in Navier-Stokes dynamics not only the velocity field and derivatives thereof are 
independently scaling fields as analyzed in ref. ||21|| , but that there are already two indepen- 



dently scaling velocity fields v L (r) and v T (r). Indeed, we will find significant differences for 
the longitudinal and transversal scaling corrections, namely, to put the result in a nutshell, 
8^ = 0.21 ± 0.01 and 6$ = 0.43 ± 0.01 in very good agreement with Boratav's result for 



decaying turbulence |TT[ 



The question which immediately comes up is whether this finding originates in the 
anisotropy of the flow. Therefore, we carefully analyze the degree of anisotropy of the 
numerical flow. We find only limited anisotropy and only on the very large scale and there- 
fore consider the different scaling of longitudinal and transversal structure functions the 
more remarkable. 

The second point we examine is whether these different intermittency corrections 8£p 
and 5£p correspond to different hierarchies of the moments. Such hierarchies were suggested 
by She and Leveque 0] for the r- averaged energy dissipation rate e r , namely 

( 6 r +l ) _ p>> ( ( e r) V / (oohl-/3 3 / fi \ 

B' p constant, = lim p ^ 00 ((e^ +1 )/(e^)); the SL parameter (3 is therefore called hierarchy 
parameter. Such a hierarchy means that the corresponding probability distribution function 
obeys a log-Poisson statistics . 

Ruiz Chavarria et al. |2]| extended the idea of hierarchies to structure functions. As- 
suming Kolmogorov's refined similarity hypothesis D p (r) ~ (e p / 3 )r p ^ 3 p, 24J the structure 



function hierarchy can be derived from the SL hierarchy (|g) [|18| , [231 and reads 

D p+ i(r) _ 



D p (r) 




(7) 



B; constant, D (oo) (r) = lim p ^ 00 ( J D p+1 (r)/ J D p (r)) = (re^)) 1 / 3 . 

We will calculate the hierarchy parameter j3 both for the longitudinal and transversal 
structure functions, very carefully considering the systematic and statistical errors. First, we 
find significant deviations between [3 L and f3 7 '. Second, we find a slight but also significant 
dependence of the hierarchy parameters (3 L and /3 T on the order of the moment which is not 
expected within the SL model. 

The paper is organized as follows: In section 2, we define the numerical flow and carefully 
check its isotropy, in section 3 we report on various scaling relations, employing ESS and 
the generalized ESS (GESS, [|l5|,[l^]); we also calculate the hierarchy parameters (3 L and j3 T 



In section 4 we determine 8£ and 8^ within a reduced wave vector set approximation of 



the Navier-Stokes dynamics [25|-g7j i n which very large Re\ can be achieved. Conclusions 
are drawn in section 5. 
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FIG. 1. Second order transversal structure functions Dj(r) for the isotropic N = 96, v = 0.006 
simulation (Re\ = 110, upper), the isotropic N = 60, u = 0.009 simulation (Re\ = 70, middle), 
and the anisotropic N = 60, v = 0.009 simulation {Re\ = 70, bottom). The solid lines are 
calculated from the definition (§) of Dj(r) (for the two directions being perpendicular to r), the 
dashed line is calculated from relation (Q) which holds for perfect isotropy and homogeneity. For 
the anisotropic case anisotropy can be seen on all scales; also the ISR slope deviates from the 
expected value (2 = 0.70. 



II. SET UP OF THE FLOW AND CHECK OF ITS ISOTROPY 

The 3D incompressible Navier-Stokes equations are numerically solved on a N 3 grid 
with periodic boundary conditions. Spherical truncation is used to reduce aliasing. For 
the isotropic flow simulation (denoted by "I") we force the system on the largest scale 
(wavevectors k = (0,0, ±1)/L, k = (0, ±1, ±1)/L, k = (±1, ±1, ±1)/L, and permutations 
thereof) with a forcing term as e.g. described in ref. |[27|| . Units are fixed by picking the 
length scale L = 1 and the average energy input rate (= the energy dissipation rate) e = 1. 
The Taylor-Reynolds number is denned as Re\ = ui^ ms \/i , ) where A = ui^ rms / (diUi) rms is 
the Taylor length and v the viscosity. Most of our results refer to iV = 96 and v = 0.006, 
corresponding to a resolution of scales r > 2ttL/N w 3i] and Re\ = 110. Time integrations 
of about 140 large eddy turnover times are performed. Averages are taken over space 
and time. To check the Reynolds number dependence we also did an isotropic N = 60, 
v = 0.009 simulation (240 large eddy turnovers) which has Re\ = 70. For a less isotropic 
flow simulation (denoted by "A") we only force one mode k = (0,0, 1)/L. This simulation 
is done for N = 60, v = 0.009, for about 210 large eddy turnovers; it has Re\ = 70, too. 

We checked the isotropy of the flow in several ways: 

• We calculated the structure functions for different space directions and compared them 
among each other. For the simulation "I" good agreement is found, for "A" one space 
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FIG. 2. Third order structure function D^(r), directly calculated from the numerics 
(long dashed) and from Kolmogorov's structure equation (||) (solid). Also shown are D^ L {r) 
(dashed-dotted) and l-D^r)! (short dashed). The data are for the isotropic N = 96 simulation 
{Re x = 110). 



direction is distinguished as expected from the type of forcing, see figure |l|. Moreover, 
for the isotropic simulation we find less than 5% deviations between (uf), and 
(«§}. Note that D\ (r = n) and Dj{r = 7r) (We used r = tc as the largest possible 
space separation in numerical flow with periodic boundary conditions.) do not equal 
2(w^) (z = 1,2, or 3) as expected for experimental, isotropic flow at r = oo. We 
find deviations up to 25% which means that the velocities are still correlated at the 
space distance of r = n. For the longitudinal velocities we find a positive correlation 
of about 25%, for the transversal velocities we find a negative correlation of about 
15%. Geometrically, this means that there is a large scale eddy with diameter ~ n. 
- We can not fully exclude that the results on scaling exponents we will report on 
are influenced by the flow geometry (periodic boundary conditions). They might be 
different for different geometries (e.g., those in experimental flows). 

• We checked relation (|4]) which only holds for isotropy |l[. For "I" there are only large 
scale deviations, for "A" deviations show up down to small scales, see figure |1|. 

• We checked the relation @, see figure 0. It holds for isotropic flow. The agreement is 
reasonable. However, there still is no developed inertial subrange due to the low Rex- 
The curve looks very similar to the experimental curve for comparable Rex, cf. fig. 
2 of ref. |§. In particular, also the experimental curves bend down for large r. The 
reason for this of course is that at large scales the fluctuations are Gaussian and odd 
order moments vanish. We ascribe the deviations in the viscous subrange (VSR) to 
the lack of perfect convergence of odd moments. This difference remained even for as 
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FIG. 3. Isotropy coefficient I(ki) eq. (|9|) for the simulations "I" (solid) and "A" (dashed) 
(N = 60). 

long averaging times as 140 large eddy turnovers. Also the relation Dj(r) = is not 
yet fulfilled for this low Re\, though the modulus of Dj (r) is more than one decade 
smaller than the modulus of D^(r) for all scales, see figure |2|. 

• For perfect isotropy, the mean energy dissipation rate e can be calculated from any 
component of the strain tensor ditij, e.g., 

1 ^ 

e = 15u((d lUl ) 2 ) = -v((d 2Ul ) 2 ). (8) 

For the isotropic flow, these relations hold very well, see table |I[ for the anisotropic 
one there are deviations up to 15%. 

• For an isotropic flow, the isotropy coefficient [EHJ 

jn E 11 (k 1 )-k 1 dE 11 (k 1 )/dk 1 fn . 

1{kl) ~ m^) (9) 

which compares the longitudinal and transversal energy spectra Eu(ki) and ^22(^1) 
should become 1 P5| . I{ki) is shown in figure |^. Indeed, in general I{ki) is closer to 
1 for "I" than it is for "A". We do not quite understand the bump around k\ = 2.5 
in I{k\) in the isotropic simulation. We tend to ascribe it to the forcing of the modes 
(±1, ±1, ±1)/L. The wiggles for very large k\ are numerical artifacts because of the 
derivative in eq. @. 
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energy input 


e 


15v((d lUl ) 2 ) 




I, N = 96 


1 


1.003 


0.984 


1.021 


I, N = 60 


1 


1.003 


0.984 


0.994 


A, N = 60 


1 


1.004 


0.931 


0.851 



TABLE II. Energy input and dissipation rates for the three numerical simulations. The good 
agreement between the energy input and the total energy dissipation rate e means statistical sta- 
tionarity. The degree of agreement between the last two columns with 1 characterizes the degree of 
isotropy in the VSR. 



III. SCALING EXPONENTS FOR LONGITUDINAL AND TRANSVERSAL 

STRUCTURE FUNCTIONS 



To determine the degree of intermittency in the longitudinal and transversal structure 
functions we employ a type of ESS HI4], |18|, |3~0|1 by calculating generalized structure functions 



and plotting them vs ("compensated ESS plot" |||31j]), see figure |[ The intermittency 
exponents <5£ 6 (the ISR slopes in figure [|) for the longitudinal and transversal structure 
functions are clearly different; the transversal signal shows considerably more intermittency. 
No dependence on Re\ is found. The values of 5£ 6 ' for the isotropic Re\ = 110 and 
Re\ = 70 simulations are given in table |. Surprisingly, also the anisotropic simulation 
"A" approximately has the same scaling exponents, namely ££g = —0.23 ± 0.01 and <5£g = 
—0.40 ± 0.01, see figure f|. Therefore, in what follows we will only focus on the isotropic 
simulation "I" . Our results for various 5£ p determined as in figure [| are summarized in figure 

It can be seen that the intermittency corrections 5£ p and 5£ p clearly deviate throughout, 
i.e., transversal velocity fluctuation are much more intermittent than longitudinal ones. 
Though it has been known for many years that the transversal velocity gradient d%Ux is 



more intermittent than the longitudinal one d\U\, see e.g. [^2] - in our simulations we have 



flatnesses of Fq 1u1 = 4.9, Fq 2U1 = 7.0, Fq iU2 _q 2Ui = 7.0 - it is not trivial that this difference, 



probing the VSR, is carried on into the ISR. In ref. [|Tl[ the different degrees of intermittency 
were associated with different types of structures: Longitudinal fluctuations with strain like 
structures, transversal fluctuations with vorticity like structures - which both is in keeping 
with the definitions of strain and vorticity, respectively. 

Let us discuss our results on 5£ p ,T in figure || in more detail. The values for 5£ p are 
well described by the SL model fit eq. (0) with the SL values C = 2, (3 = (2/3) 1 / 3 w 0.874 
as found for many other isotropic, even low Re\ number, experimental or numerical flows 
|2|, [T3|, [TJ| . For the physical interpretation of the parameters in the SL model we refer to refs. 



|T^.|22|. From a phenomenological point of view, one could consider eq. (||) simply as a two 
parameter fit of the £ p 's. The two SL parameters for the transversal scaling exponents can 
be viewed as a simple way to quantify the degree of intermittency. 
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FIG. 4. Compensated ESS type plots for D%/(D% L ) 2 vs Df (circles) and D%/(Df) 2 vs Df 
(stars). The ISR slopes are respectively (5£g and 5t^. The upper part of the figure refers to the 
isotropic simulation with N = 96, the lower one to the anisotropic one with N = 60. 
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96, squares N = 60) and 6^z (crosses 



FIG. 5. 

N = 96, triangles = 60) from the isotropic numerical simulations. The dashed lines are 1-pa- 
rameter fits within the SL-model where the /?'s have been taken fixed from the averaged fits of 
the GESS type plot, i.e., (5 L = 0.947 and f3 T = 0.870. The one free fit parameter is thus C . We 
obtain the shown remarkably good fits with Cq = 9.3 for the longitudinal data (short dashed) and 
Cj = 3.7 for the transversal data (long dashed). The standard SL fit Q = (2/3) 1 / 3 , C = 2 
fits the longitudinal corrections also pretty well, see the solid line. 



We now suggest a method to replace this one two-parameter-fit by two one-parameter 
fits. To do so, we employ generalized extended self similarity (GESS, flT5] , |I%|l ) and plot G p 
vs G q , see figure The slope p VA of such a plot is by definition fllCf ) 

For fixed p, q, p£ and p^ q are significantly different, see table |ITT[ We checked this result 
very carefully. The small error bars in p VA result from linear regressions in GESS plots as 
in figure ^ and in addition from an averaging over the different space directions. We also 
checked this for smaller Re\ = 70; the deviations in comparison to the results given in table 



[111 are never larger than 0.5%. Moreover, to make sure that our numbers are well converged, 
we also averaged over only 30, 60, and 90 large eddy turnovers rather than 140; still, the 
result is the same; the deviations are smaller than the error bars. 

To quantify the quality to which GESS holds we checked the relation p ps = p p , q p qtS , 



implied by GESS, for various p, q. Table |T| allows the reader to do so. Neither for the p^ 
nor for the did we find a single example where there were deviations larger than the 
error bars. E.g., p\ 4 ■ p\ 6 = 0.1177 ± 0.0005 which equals p\ 6 within the error bars. 

The error bars up to now stem from statistics. One would like to be able to judge the 
size of the systematic errors. Therefore, in figure [7| we display the local slope of the curves 
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FIG. 6. GESS type plot G^ 6 \r) vs G^ 2 '(r) for both the longitudinal (solid) and the transver- 
sal (dashed) G-structure functions for the isotropic simulations. The slopes of these curves are 
Pq 2 = —8.49 ±0.02 and p§ 2 = —7.52 ±0.04, respectively. The errors result from a linear regression 
of every single curve, from weighted averaging of the results for different space directions, and 
different Re\. The hardly distinguishable lines within the two bunches of curves are the results 
for different directions and different Reynolds numbers Re\ = 110 and Re\ = 70. The good agree- 
ment within the bunches means good isotropy and independence of the scaling exponents from 
Re\. Note that in this type of plot the far VSR collapses into the upper left point of these curves. 
The filled bullets refer to r = 10rj for Re\ = 110 (left one) and Re\ = 70 (right one). The open 
bullet refers to the outer length scale r = L. 



in fig. |5|. Both the longitudinal and the transversal local slope slightly increase (modulus- 
wise) with increasing scale (i.e., from left to right), which shows the limitations of the above 
statement that GESS is fulfilled with remarkable quality. The error bar calculated from 
averaging the local slope is much bigger than above statistical error. From averaging up 
to the scale r = 2.0 « IOO77 we obtain p£ 2 = -8.36 ± 0.14 and p£ 2 = -7.36 ± 0.26. The 
numbers for the p^'J from table |TTT| are within the (now about ten times larger) error bars. 
These systematic errors are one order of magnitude bigger than the purely statistical ones 
in table [[II]. We summarize the values of p^'J (and their error bars) determined in this way 
in table [IV]. Note, however, that the deviations between and pj and correspondingly 
also between the resulting /3's (see below) are still statistically significant. 

Within the SL model, the p p , g 's only depend on (3, not on any other parameter, 

= (i-/y)-(p/3)(i-i 3 ) 



For each p VA we calculate @ Piq , resulting from equation (0), and its error, see table [TTT] and 



table |TV|- If the SL model were exact, (3 should not depend on p and q. 



11 



-7 




-0.22 -0.18 -0.14 

ig G (2) 

FIG. 7. Local slopes of figure |6| for the Re\ = 110 simulation. The two dashed lines are for 
the two different space direction for the transversal structure functions, the three solid lines are for 
the three space directions of the longitudinal structure functions. If we calculate the average (for 
scales up to r = 2.0 ~ IOO77) we obtain pg 2 = —8.36 ± 0.14 and p^ 2 = —7.36 ± 0.26. The arrows 
refer to IO77 and L, respectively. 



In figure |8| we offer a 3D plot of (3^ together with the error bars resulting from the 
(larger) systematic errors of p, cf. table [IV]. From figure § the difference between (3 L and 
(3 T seems to be significant. This result is at variance with Camussi and Benzi's [^UJ who 
obtained that the difference of both (3 L and (3 T to the SL value (3 = (2/3) 1/3 « 0.874 is at 
most 1.2% = 0.010. 

Another feature of figure § is that (3^ shows a small trend towards smaller values for 
larger p, q which is not expected within the SL model. If we average over all f3 p>q nevertheless, 
we obtain f3 L = 0.947 and p T = 0.870. 

Knowing f3, there is only the parameter Cq left in eq. (^). If we take the above mean 
values (3 L = 0.947 and (3 T = 0.870, we obtain as best fits to the £ p data in figure | Cq = 9.3 
(the x 2 °f the fit is x 2 — 10) and Cq = 3.7 (with x 2 — 1)> excellently describing the 
numerical data. We do not ascribe any physical meaning to the parameter values obtained 
in our fit. Note that for our £ L data this fit is superior to the SL model with the original 
parameter values (3 = (2/3) 1 / 3 , C = 2. If we choose the SL- value (3 L = (3 T = (2/3) V 3 we 
obtain C L = 1.97 with X 2 = 10 3 and C? = 3.9 with X 2 = 0.8. (The ^-values for C L are 
larger than those for Cq as the errors of the 5£r are smaller than those of 5£T.) 

We now directly check the hierarchies of the structure functions [12|, [18], |23| . From eq. 
(0) it is easy to derive [P~8"|,|2"3" 



F p+1 (r) = B p (F p (r)fF(r) (13) 

with 
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P,Q 


L 
Pp,q 


qL 
p 


T 
Pp,q 




2,4 


-0.5093 ± 0.0004 


0.973 ± 0.001 


—0.5440 ± 0.0011 


0.877 ± 0.003 


2,5 


-0.2073 ± 0.0002 


0.964 ± 0.001 


—§?IT12 ± 0.0006 


0.875 ± 0.002 


2,6 


— 0.H76 ± 0.0002 


0.957 ± 0.001 


— 0.1318 ± 0.0004 


0.873 ± 0.002 


2,7 


-0.0774 ± 0.0001 


0.952 ± 0.001 


-0.0884 ± 0.0003 


0.872 ± 0.002 


2,8 


-0.0555 ± 0.0001 


0.947 ± 0.001 


-0.0645 ± 0.0003 


0.871 ± 0.003 


4,5 


0.4072 ± 0.0001 


0.947 ± 0.001 


0.4177 ± 0.0003 


0.870 ± 0.002 


4,6 


0.2311 ± 0.0001 


0.941 ± 0.001 


0.2423 ± 0.0004 


0.869 ± 0.003 


4,7 


0.1521 ± 0.0002 


/ \ s \ i \ t—? i r\ r\r\-\ 

0.937 ± 0.001 


0.1625 ± 0.0005 


0.868 ± 0.003 


4,8 


0.1092 ±0.0003 


0.933 ± 0.002 


0.1186 ±0.0005 


0.867 ± 0.003 


5,6 


0.5675 ± 0.0002 


0.935 ± 0.001 


0.5800 ± 0.0005 


0.867 ± 0.003 


5,7 


0.3736 ± 0.0004 


0.931 ± 0.002 


0.3890 ± 0.0008 


0.867 ± 0.003 


5,8 


0.2682 ± 0.0006 


0.928 ± 0.002 


0.2840 ± 0.0010 


0.866 ± 0.004 


6,7 


0.6583 ± 0.0005 


0.927 ± 0.002 


0.6708 ± 0.0008 


0.865 ± 0.004 


6,8 


0.4725 ± 0.0008 


0.924 ± 0.003 


0.4901 ± 0.0014 


0.863 ± 0.005 


7,8 


0.7177 ±0.0007 


0.922 ± 0.003 


0.7308 ± 0.0010 


0.860 ± 0.005 



TABLE III. p^'J for various pairs p, q from GESS type plots as in figure 6. By definition, 
Pp,q — Pq,p- The errors are purely statistical ones. In the third and fifth column, we give the j3's 
resulting from eq. (12). Note that GESS implies p PjS = p p ,qPq )S - This relation can be used to check 
the quality of GESS. 



_ D p+1 (r) 

Fp+1 ~ (14) 

and 

/ n (r) x^-v/w-p 3 )] 

p{ ^(im^) • (15) 

With the mean /3's obtained above, we plot F p+1 (r) vs (F p (r)) /3 F(r). If eq. ([13]) and equiv- 
alently eq. (0) hold, the slope should be 1. Indeed, the slope is very close to 1, see figure |S| 
and table 0, which gives further support for (3 L and (3 T being different. The best agreement 
is found for p = 4 to p = 6, see table [V|. The reason is that the mean /3's best agree with 
(3 P)<1 if p, q are around 4-6, see table |ITT|. For the other p one could improve the fit by using 
the corresponding /3 p>g ; however, note that the sixth order structure function always enters 
via F(r), cf. eq. (|1~5|). 



Moreover, we find a p dependence of the prefactor B p in eq. (|13l) . Therefore, determining 
j3 from eq. ( |T3"D by plotting lgF p+1 (r) vs lgF p (r) for fixed r as a function of p as done in 
refs. |||33|] does not seem to be possible here. 
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P,Q 


L 
Pp,q 


p 


T 
Pp,q 


pT 
P 


2,4 


r\ t* -\ A i r\ r\r\ t* 

—0.514 ± 0.005 


0.958 ± 0.016 


—0.556 ± 0.012 


r~\ n < a i a aoi 

0.846 ± 0.031 


2,5 


—0.210 ± 0.003 


0.950 ± 0.014 


-0.234 ± 0.007 


0.848 ± 0.027 


2,6 


c\ i c\r\ i r\ r\r\c\ 

—0.120 ± 0.002 


0.944 ± 0.012 


—0.136 ± 0.005 


a o /in i n nor 

0.849 ± 0.025 


2,7 


r\ r\<—rr\ i r\ r\r\c\ 

—0.079 ± 0.002 


0.940 ± 0.012 


—0.092 ± 0.003 


0.851 ± 0.023 


2,8 


—0.057 ± 0.001 


0.936 ± O.Oll 


—0.067 ± 0.003 


0.852 ± 0.021 


4,5 


0.409 ± 0.001 


0.936 ± O.Oll 


0.420 ± 0.003 


0.851 ± 0.021 


4,6 


0.233 ± 0.002 


0.931 ± 0.010 


0.245 ± 0.003 


r\ n r o i r\ c\c\r\ 

0.853 ± 0.020 


4,7 


0.153 ± 0.002 


0.927 ± 0.010 


0.165 ± 0.003 


0.855 ± 0.018 


4,8 


0.110 ±0.001 


0.924 ± O.Oll 


0.120 ±0.003 


0.856 ± 0.018 


5,6 


0.569 ± 0.002 


0.926 ± 0.010 


0.582 ± 0.003 


0.855 ± 0.018 


5,7 


0.376 ± 0.002 


0.923 ± 0.010 


0.391 ± 0.004 


0.857 ±0.018 


5,8 


0.270 ± 0.003 


0.920 ± O.Oll 


0.286 ± 0.004 


0.858 ± 0.018 


6,7 


0.660 ± 0.002 


0.9l9± O.Oll 


0.672 ± 0.003 


0.858 ± 0.018 


6,8 


0.475 ± 0.004 


0.917 ±0.012 


0.491 ± 0.005 


0.859 ± 0.018 


7,8 


0.719 ±0.003 


0.914 ±0.014 


0.731 ± 0.004 


0.860 ± 0.019 



TABLE IV. for various pairs p, q determined from the local slope of the GESS type plots, cf. 
figure 7. The errors are the systematic ones, stemming from the local slope not being constant. In 
the third and fifth column, we again give the (3's and their errors resulting from eq. (12). Averaging 
the (3 's determined this way gives 

(3 L = 0.930 and (3 T = 0.855. 



IV. SCALING RELATIONS WITHIN REWA 



Very large Re\ in numerical turbulent flow can be achieved in the reduced wave vector 



set approximation (REWA) of the Navier-Stokes equation [[25H 271, |34], |35| . REWA uses a 



reduced, geometrically scaling subset of wavevectors on which the Navier-Stokes equation is 
solved. Here we choose a basic set of 50 modes per level. Very high Taylor- Reynolds numbers 
up to Re\ = 7 ■ 10 4 [^7|,|34[] can be achieved, however, flow structures are underrepresented 
31| and the intermittency corrections are strongly underestimated [pT , |35[ . 

We redid ESS types plots for REWA for Re x = 8 • 10 2 and for Re x = 1.4 • 10 5 for both 
the longitudinal and the transversal sixth order structure functions, see figure ID. There is 
no detectable difference between the longitudinal and transversal scaling exponents. The 
absolute value 5^ ~ 5£j ~ —0.009 is much smaller (modulo wise) than the experimental 
or above numerical value <5£g ~ —0.21, as extensively analyzed and discussed in the pre- 
vious work on REWA p5|-p7|, |34|, |35H . Note however that the relative error for the ££'s is 
much larger than in the full numerical simulations - we cannot exclude different degrees of 
intermittency for and within REWA. 

We do not know whether our results on REWA indicate that the differences between 
the scaling of Dp vs D% T and vs D% L observed in the above full numerical simulations 
for small Re\ are finite Re\ effects or whether they are artefacts of the REWA thinning of 
large wave vectors [^5],^7], connected to the suppression of small scale structures. Such 
structures are associated with the different scaling of longitudinal and transversal structure 
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0.95 




FIG. 8. The longitudinal and transversal hierarchy parameters Pp q (upper) and f3 p r „ (lower) 



respectively. The data are taken from table III, i.e., for the "I" simulation with Re\ = 110, the 



averaging time is 140 large eddy turnovers. To get an idea of the size of the error, two error bars 
are drawn, representing the much larger systematic errors rather than the statistical ones. First, 
we observe that (3 L and (3 T are clearly different. Second, a slight dependence of (3 L ' T on p, q is 
seen. 



functions in ref. []n}| . As within REWA no Re\ dependence of the 5£ and 5£ is observed, 
see figure [K], we favor the second interpretation. 

We also tried GESS type scaling within REWA. No statistically significant deviations 
between the p^ q and the p^ q were found. 



V. SUMMARY, CONCLUSIONS, AND OUTLOOK 

To summarize, we offer strong evidence that the transversal velocity fluctuations show 
stronger intermittency than the longitudinal ones. Our numerical values for the longitudinal 
and the transversal scaling exponents £^ and £j for forced stationary turbulence agree very 
well with those of Boratav and Pelz for decaying turbulence Jll| , see table |. This finding is 



independent of Re\, at least for the relatively low Re\ we examined. For an anisotropic flow 
we essentially obtained the same scaling exponents. Only for the REWA calculations which 
underrepresent the small scale structures of the flow we do not find a statistically significant 
deviation between 8^ and <5£j, however, the relative error is much bigger than for the full 
simulations. 

We reiterate that Q = (2 because of relation ([|); a generalization of this equality to 
higher order moments p > 2 is wrong. 

GESS is fulfilled with satisfactory precision for both longitudinal and transversal struc- 
ture functions. The GESS scaling exponents p^ q and p^ q are different. This result is the 
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FIG. 9. Log-Log plot of F p+ i(r) vs (F p (r))^F(r) for various p for the longitudinal and transver- 
sal structure functions, (3 L = 0.947, f3 T = 0.870. The lines are arbitrarily shifted in order for the 
slopes to be visible. The upper 5 lines are for the longitudinal structure functions, p = 3 (upper) 
to p = 7 (lower), the lower 5 lines are for the transversal structure functions, p = 3 (upper) to 
p = 7 (lower). 



more remarkable, as those of the longitudinal velocity structure function agree with those 
of an active [BH] or passive scalar |T3 or even with those calculated for the magnetic field in 



MHD [0, see table 1 of ref. @ 



The She-Leveque hierarchy parameters j3 L)T following from the //^-exponents conse- 
quently also differ. Both show a weak dependence on p, q which is not expected within the 
She-Leveque model. 

In the whole analysis we took great care of systematic and statistical errors to get sig- 
nificant statements. 

To conclude, there seem to exist independently scaling velocity fields v L (r) and v T (r), 
i.e., the Navier-Stokes dynamics seems to make use of this degree of freedom being allowed 
by symmetry. It is very likely that the two different scaling velocities fields will also be 
reflected in the flow geometry. 

A more complete discription of the statistics of the velocity field was recently suggested 
by L'vov, Podivilov, and Procaccia ||38|| . These authors point out that because of rotational 
symmetry only the SO (3) irreducible amplitudes of the velocity structure tensors should 
obey clean scaling. In ref. we analyse the scaling properties of these amplitudes numer- 
ically and show (for fourth order moments) how they are connected to the longitudinal and 
transversal structure functions. 

Presently, neither the longitudinal nor the transversal scaling exponents, nor the scaling 
exponents of the irreducible SO (3) invariants of the velocity correlations [Q can be calcu- 
lated analytically from the Navier-Stokes dynamics. Many phenomenological models based 
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p 


/3 L = 0.947 
slope L 


B p 


ff = 0.870 
slope T 


B p 


f3 L = (2/3) 1 / 3 
slope L 




3 


0.9984 ± 0.0001 




0.975 


0.9988 ± 0.0002 


0.981 


0.9988 ± 0.0001 


0.981 


4 


0.9995 ± 0.0001 




0.944 


0.9993 ± 0.0004 


0.961 


0.9936 ± 0.0001 


0.969 


5 


1.0020 ± 0.0004 




0.925 


1.0000 ±0.0005 


0.951 


0.9895 ± 0.0002 


0.967 


6 


1.0050 ± 0.0007 




0.910 


1.0020 ±0.0007 


0.945 


0.9854 ± 0.0004 


0.970 


7 


1.0080 ± 0.0010 




0.899 


1.0040 ±0.0016 


0.940 


0.9810 ± 0.0007 


0.977 



TABLE V. Slopes o/lgF p+1 (r) vs\g((F p (r))^F(r)), varying r. According to eq. (13), the slopes 
should be 1, which is pretty well fulfilled for [3 = 0.947 and (3 T = 0.870. For comparision, we also 
give the slopes if the SL-value f3 L = (2/3) ^ is used for the longitudinal structure function. The 
deviations of the slope to 1 are larger. The constants B p in (13) show a slight p- dependence. 

on various views on how intermittency develops are able to fit the longitudinal intermittency 
corrections. It should be possible to derive also the transversal intermittency corrections in 
the framework of the thinking these models are based on in order to check their consistency. 
The ultimate goal, however, must be to derive both longitudinal and transversal scaling 
exponents from the Navier-Stokes equation. 
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